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T— I Abstract 

o 

Human motion in the vicinity of a wireless link causes variations in the link received signal strength 

use these RSS variations in a static wireless network to detect, locate and track people in the area of 

the network, even through walls. However, intrinsic motion, such as branches moving in the wind and 
I p rotating or vibrating machinery, also causes RSS variations which degrade the performance of a DFL 

system. In this paper, we propose and evaluate two estimators to reduce the impact of the variations 
Q caused by intrinsic motion. One estimator uses subspace decomposition, and the other estimator uses a 

least squares formulation. Experimental results show that both estimators reduce localization root mean 
" squared error by about 40% compared to VRTI. In addition, the Kalman filter tracking results from both 

estimators have 97% of errors less than 1.3 m, more than 60% improvement compared to tracking results 

from VRTI. 

d 

, I I. Introduction 

As an emerging technology, device-free localization (DFL) using radio frequency (RF) sensor networks 
^ has potential application in detecting intruders in industrial facilities, and helping police and firefighters 

track people inside a building during an emergency lUl. In these scenarios, people to be located cannot be 
expected to participate in the localization system by carrying radio devices, thus standard radio locahzation 
techniques are not useful for these applications. 

Various RF measurements including ultra-wideband (UWB) and received signal strength (RSS) have 
been proposed and applied to detect, locate and track people who do not carry radio devices in an indoor 
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environment ||2l, ||3l, H, Q. RSS measurements are inexpensive and available in standard wireless 
devices, and have been used in different DFL studies with surprising accuracy H, ||6l, 151. These RSS- 
based DFL methods essentially use a windowed variance of RSS measured on static links. For example, 
|[5l deploys an RF sensor network around a residential house and uses sample variance during a short 
window to track people walking inside the house; ||6l places RF sensors on the ceiling of a room, and 
track people using the RSSI dynamic, which is essentially the variance of RSS measurements, with and 
without people moving inside the room. In this paper we use windowed variance to describe the various 
functions of RSS measurements recently used in different DFL studies H, Q, Q, 13, and we call these 
methods variance-based DFL methods. 

For variance-based DFL methods, variance can be caused by two types of motion: extrinsic motion 
and intrinsic motion. Extrinsic motion is defined as the motion of people and other objects that enter 
and leave the environment. Intrinsic motion is defined as the motion of objects that are intrinsic parts 
of the environment, objects which cannot be removed without fundamentally altering the environment. 
If a significant amount of windowed variance is caused by intrinsic motion, then it may be difficult to 
detect extrinsic motion. For example, rotating fans, leaves and branches swaying in wind, and moving or 
rotating machines in a factory all may impact the RSS measured on static links. Also, if RF sensors are 
vibrating or swaying in the wind, their RSS measurements change as a result. Even if the receiver moves 
by only a fraction of its wavelength, the RSS may vary by several orders of magnitude. We call variance 
caused by intrinsic motion and extrinsic motion, the intrinsic signal and extrinsic signal, respectively. 
We consider the intrinsic signal to be "noise" because it does not relate to extrinsic motion which we 
wish to detect and track. 

This work is motivated by our inability to achieve the performance of 0.6 m average tracking error 
reported in [5 1 in a repeat of the identical experiment in May, 2010. Our new experiment was performed 
at the same location and using the identical hardware, number of nodes, and software. Yet, in the new 
experiment, variance-based radio tomographic imaging (VRTI) does not always locate the person walking 
inside the house as accurately as reported in O. Sometimes the position estimate error is as large as 
six meters, as shown in Figure |6] Investigation of the experimental data quickly indicates the reason for 
the degradation: periods of high wind. Consider the RSS measurements recorded during the calibration 
period, when no people are present inside the house. From the calibration measurements of iQ, the 
standard deviations of RSS measurements are generally less than 2 dB. However, the RSS measurements 
from our May 2010 experiment are quite variable, as shown in Figure [T] The RSS standard deviation 
can be up to 6 dB in a short time window. Considering there is no person moving inside the house, that 
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is, no extrinsic motion during the calibration period, the high variations of RSS measurements must be 
caused by intrinsic motion, in this case, wind-induced motion. 
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Fig. 1: Intrinsic signal measurements: RSS measurements from three links during the calibration period 
when no people are present in the environment. 



The variance caused by intrinsic motion can affect both model-based DFL and fingerprint-based DFL 
methods. To apply various DFL methods in practical applications, the intrinsic signal needs to be identified 
and removed or reduced. The subspace decomposition method has been used in spectral estimation, sensor 
array processing, and network anomaly detection lUl, |l9l, lITOll . iFTTl . We apply this method to VRTI, 
which leads to a new estimator we refer to as subspace variance-based radio tomography (SubVRT) ||T2|. 
Inspired by the fact that SubVRT makes use of the covariance matrix of link measurement and significantly 
reduces the impact of intrinsic motion, in this paper, we formulate a least squares (LS) solution |]T3l for 
VRTI which uses the inverse of the covariance matrix. We call this method least squares variance-based 
radio tomography (LSVRT). 

The contribution of this paper is to propose and compare two estimators - SubVRT and LSVRT to 
reduce the impact of intrinsic motion in DFL systems. Experimental results show that both estimators 
reduce the root mean squared error (RMSE) of the location estimate by more than 40% compared to 
VRTI. Further, we use the Kalman filter to track people using localization estimates from SubVRT and 
LSVRT. The cumulative distribution functions (CDFs) of the tracking errors show that the tracking results 
from SubVRT have 97% of errors less than 1.4 m, a 65% improvement compared to VRTI, while 97% 
of tracking errors from LSVRT are less than 1.2 m, a 70% improvement. 

The rest of this paper is organized as follows: Section [IT] discusses the subspace decomposition method 



and least squares method for noise reduction in DFL. Section III describes the experiments. Section IV 
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shows the experimental results, and Section |V] investigates the Kalman filter tracking. Related work is 



presented in Section NT\ and the conclusion is given in Section VII 



II. Methods 

In this section, we formulate a variance-based DFL problem, introduce the subspace decomposition 
method, and propose our SubVRT estimator. After that, we use the measurement covariance matrix in a 
least squares (LS) formulation and propose another estimator, LSVRT. Finally, we discuss the connection 
between these two estimators. 



A. Problem statement 

For an RF sensor network with N sensors (radio transceivers) deployed at static locations, we use 
to denote the coordinate of sensor j. Each sensor makes an RSS measurement with many other sensors, 
and we use s/ ^ to denote the RSS measured at node ii sent by node ji at time t, where ii and ji are 
the receiver and transmitter number for link /, respectively. Time t is discretized, thus t G Z. We assume 
constant transmitter power so that changes in ^ are due to the channel, not to the transmitter. Then we 
denote the windowed RSS variance as: 

^ m— 1 

yi,t = 7 {si^t - si^t-if (1) 

i=0 

where m is the length of the window, and j = ^ X]2=o^ si,t-i is the sample average in this window 
period. 

Consider that the network has L directional links on which we measure signal strength (in general, 
L < N{N — 1)). We let y^*) = [yi^t,y2,t, ■ ■ ■ , yL,t\^ be the vector of windowed RSS variance from all 
L links at time t. If we do not need to represent time, we simplify the notation to y = [yi, 2/2, • " " ) yiY ■ 
Then we use yc to denote the calibration measurements collected during the calibration period, when no 
people are present in the environment; and we use y^ to denote the measurements from the real-time 
operation period. The goal of DFL is to locate people during real-time operation. 

For VRTI, a model-based DFL method, the presence of human motion within P voxels of a physical 
space is denoted by x = [xi, X2, xp]"^, where x, = 1 if extrinsic motion occurs in voxel i, and Xj = 
otherwise. Work in [5J has shown the efficacy of a linear model that relates the motion image x to the 
RSS variance jr'- 

y^ = VFx + n (2) 
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where n is an L x 1 noise vector including intrinsic motion and measurement noise, and W is an L x P 
matrix representing the weighting of motion in each voxel on each link measurement. The weighting of 
voxel p on link / is formulated as llH : 

1 di,^p + dj,^p < di,j, + X 

Wi^p = -j== I (3) 

V^»'J' I otherwise 

where is the Euclidean distance between two sensors ii, ji on link / located at and j, ; dj^^p is 
the Euclidean distance between sensor ji and Zp, the center coordinate of voxel p; di^.p is the Euclidean 
distance between sensor ii and voxel p; A is a tunable parameter controlling the ellipse width, and is 
a constant scaling factor. 

Once we have the forward model, the localization problem becomes an inverse problem: to estimate 
P dimensional position vector x from L dimensional link measurement vector yr. Certain regularization 
methods are necessary for this ill-posed inverse problem, and it is shown in Q that submeter locaUzation 
accuracy can be achieved by using the Tikhonov regularization. Thus, we use the Tikhonov regularized 
VRTI solution, which is given as: 

Hi = {W^W + aQ^Q)-^W^ (4) 
where Q is the Tikhonov matrix, and a is a regularization parameter. 



B. Subspace decomposition method 

The subspace decomposition method has been widely used in spectral estimation, sensor array process- 
ing, etc. im, lOn to improve estimation performance in noise. It is closely related to principal component 
analysis (PC A), which is widely used in finding patterns in high dimensional data fT?!. 

From the L-dimensional cahbration measurement vectors yc, we may estimate its covariance matrix 
Cye as: 

A/-1 

Cy. = ^ j:(y(*)-/.e)(y^*^-/^cf (5) 

i=0 

where M is the number of sample measurements, yi*^ is the calibration measurement vector yc at time 
t, fi^= jj Y^tio^ yc'' is the sample average. 

Then, we perform singular value decomposition (SVD) on Cy^: 

Cy^ = UAU^ (6) 
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where the unitary matrix U = [ui,--- ,ul], and the diagonal matrix A = diag {Ai, A^}. Right 
multiplying U on both sides of Q, we have: 

Cy^Ui = XiUi (7) 

where Uj is the eigenvector corresponding to the eigenvalue Aj. If the eigenvalues are in descending 
order, the first principal component ui points in the direction of the maximum variance in the calibration 
measurements, the second principal component U2 points in the direction of the maximum variance 
remaining in the measurements, and so on. If the first few eigenvalues are much larger than the others, 
then most of the variance in the calibration measurements can be captured by these principal components. 

We perform the above PCA procedures on calibration measurements from two sets of experiments as 
described in Section III The eigenvalues of Cy^ from these experiments are shown in Figure [2] Because 
there is more intrinsic motion in Experiment 2, we see that the largest eigenvalue from Experiment 2 
is almost twice as large as that from Experiment 1. We also see that for Experiment 1, the first four 
eigenvalues are much larger than the other eigenvalues, thus the corresponding eigenvectors can capture 
most of the variation in the measurements. However, for Experiment 2, there are more large-valued 
eigenvalues, and more eigenvectors are necessary to represent the major variation in the measurements. 



160 




Component number 



Fig. 2: Scree plot. 



From the scree plot, we decide how many principal components, k, are necessary to capture the 



majority of the variations (we discuss selection of k in more detail in Section IV-Cl. Then, we use 
a lower dimensional space spanned by these principal components to represent the space containing 
the majority of the intrinsic signal measurements. Thus, in subspace decomposition, we divide all the 
principal components into two sets: U = [ui,U2,-- - ,Ujt] and U = [u^+i, • " " lU^]- Since the 



October 10, 2011 



DRAFT 



7 



variance during the calibration period is caused by intrinsic motion, that is, the variance captured by U is 
intrinsic signal, we call the subspace spanned by U the intrinsic subspace, and the other subspace spanned 
by U the extrinsic subspace. Once the two subspaces are constructed, we can decompose the measurement 
vector y into two components - intrinsic signal component y and extrinsic signal component y: 

y = y + y- (8) 

Since the principal components are orthogonal, the intrinsic signal component y and the extrinsic signal 
component y can be formed by projecting y onto the intrinsic subspace and the extrinsic subspace, 
respectively: 

y = njy = UU^y (9) 
y = Hey = (/ - UU^)y (10) 

where 11/ = UU'^ is the projection matrix for the intrinsic subspace, and He = / — II/ is the projection 
matrix for the extrinsic subspace. 

The key idea of SubVRT is to use the decomposed extrinsic signal component of the measurements in 
VRTI. We project the real-time measurement vector y^ onto the extrinsic subspace to obtain the extrinsic 
signal component y-r = {I — UU'^)yr- Then, we replace y^ in ^ with y^ and obtain the solution of 
SubVRT: 

n2 = {W'^W + aQ^Q)-^W^{I -UU^). (11) 



From ( [TT) , we see that the solution is a linear transformation of the measurement vector. The trans- 
formation matrix 112 is the product of the transformation matrix Hi in (|4]) with the projection matrix 
for the extrinsic subspace II^;: 112 = nill^;. Since the transformation matrix 112 does not depend on 
instantaneous real-time measurements, it can be pre-calculated, and it is easy to implement SubVRT for 
real-time applications. Calculation of x from y^ requires LP multiplications and additions. 



C. Least squares method 

SubVRT performs SVD on the calibration measurement covariance matrix. Here, we introduce our 
LSVRT estimator formulated as a least squares (LS) solution, which uses the inverse of the covariance 
matrix. 
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1) Formulation: To derive the least squares solution to the linear model expressed in (|2]), the cost 
function can be written as Ill3l: 



J(x) = ||VFx-y^||^_^ + ||x-xa||^^ (12) 
= (y, - W^YC-\yr - W^) + (x - Xa)^C, ^(x - xj 

where ||n||^ indicates weighted quadratic distance n^Cn^n, Cn is the covariance matrix of n, Xq is 
the prior mean of x, and Cx is the covariance matrix of x. 



Taking the derivative of (12i and setting it to zero results in the LSVRT solution: 



^LS = {W^C-^W + i)-i(W^^C„ Vr + ^x,). (13) 

Since the prior information x^ can be included in the tracking period, here we assume x^ is zero, then 
([13]) becomes: 

Hs = {W^C-^W + C^^Y^W^C:^^. (14) 

The LSVRT formulation can be also justified from a Bayesian perspective. If we assume y^ conditioned 
on X is Gaussian distributed with mean H^x and covariance matrix Cn, and x is Gaussian distributed with 
mean x^ and covariance matrix Cx, then maximizing the posteriori distribution p(x|yr) is equivalent to 
minimizing the cost function in ( 12l, thus the maximum a posteriori (MAP) solution is the same as ( [T3] ). 



2) Covariance matrix Cn-" From the LSVRT solution ( [13] ), we see that the inverse of the covariance 
matrix Cn (a.k.a., the precision matrix) is needed. We may use the sample covariance matrix if the sample 
size M is greater than the number of link measurements L. However, for an RF sensor network with L 
directional links, M is typically less than L. Thus, for high dimensional problems, the sample covariance 
matrix becomes an ill-posed estimator, it cannot be inverted to compute the precision matrix. 

For high dimensional covariance matrix estimation problems, many types of regularized covariance 
matrix estimators have been proposed |[T5l . |[T6l . Here, we use the Ledoit-Wolf estimator, which is a 
linear combination of the sample covariance matrix and a scaled identity matrix, and is asymptotically 
optimal for any distribution HI 51 : 

Cn = yi^I + (1 - y)C*^ (15) 

where C* is the sample covariance matrix, ^ is the scaling parameter for the identity matrix /, and v is the 
shrinkage parameter that shrinks the sample covariance towards the scaled identity matrix. Since there is 
no extrinsic motion during calibration period, that is, x = 0, thus yc = n, and we approximate C* = Cy^. 
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Then we follow |[T5l to calculate parameters v and /i. From the Bayesian perspective, this covariance 
matrix estimator can be seen as the combination of the prior information and sample information of the 
covariance matrix. 

3) Covariance matrix C^-' The LSVRT solution also requires the covariance matrix Cx- As a means 
to generate a general statistical model for Cx, we assume that the positions of people in the environment 
can be modeled as a Poisson process. Poisson processes are commonly used for modeling the distribution 
of randomly arranged points in space. 

Analysis of Poisson point processes leads to a covariance function that is approximately exponentially 
decaying lITTl . and the exponential spatial covariance model is shown to be effective to locate people 
in an RF sensor network [18|. Thus, in this paper, we use an exponentially-decaying function as the 
covariance matrix of the human motion. 

2 /II II 



where cr^. is the variance of the human motion, 6 is a space constant, and ||xj — XjH^^ is the Euclidian 
distance between Xj and Xj. 



D. Discussion 

The SubVRT estimator and the LSVRT estimator are closely related. LSVRT needs to calculate the 
inverse of the covariance matrix Cn, while SubVRT needs to perform SVD on the sample covariance 
matrix Cy^. In this section, we show connections between these two estimators. 



First, for SubVRT, once we choose the parameter k, we can find a diagonal matrix S = diag < 



0,0,-- - ,0, ,1 



k 

such that USU^ = I — UU^. Then, the project matrix for the SubVRT solution can be rewriten as: 



Ha = {W W + aQ' Q)~^W^USU' . (17) 



For the LSVRT solution (14i and the Ledoit-Wolf covariance estimator in (15 1, if we approximate 



C* = Cy^, then the inverse of Cn can be written as 

Vfi 1 — 1/ 



C-' = —I+- -Cy!- (18) 



Substituting (j6]l in ( 18 1, we express C„ ^ in terms of A: 

= C/ci(A-i + C2/)C/^ (19) 
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where ci = and C2 = Replacing the second C^^ in (Il4|) by (19 1, the project matrix for the 



UfJ, 

LSVRT solution becomes: 



Hs = (W^C-^W + C^^)-^W'^Uci{A-^ + C2l)U^. (20) 



Now we compare the two projection matrices ( 17 1 and (|20|) in the SubVRT and LSVRT solutions. 



From the latter part of (20 1, we see that LSVRT uses ci(A ^ + C2I) to give less weights to the linear 



combinations of measurements in the eigen-space with high variance (large eigenvalues). For SubVRT, 



the diagonal matrix 5 in ( 17 1 is used to directly remove eigenvectors that correspond to the first k largest 



eigenvalues. From the former part of (17i and (20 1, we see that the inverse of the covariance matrix 



Cx ^ in the LSVRT solution plays the same role of regularization as the term aQ^Q in the SubVRT 
solution. We also see that the LSVRT estimator includes the precision matrix ^ as a weight matrix in 
W'^C^^W, while the SubVRT estimator just uses W^W. 

III. Experiments 

We use measurements from two sets of experiments in this paper. We use the data set from the 
measurements conducted in March, 2009 reported by ||5j|. We call this data set Experiment 1. The second 
experiment is a new experiment performed in May, 2010 at the same residential house, which we call 
Experiment 2. In both experiments, thirty-four TelosB nodes are deployed outside the living room of the 
house. As shown in Figure [4j eight nodes are placed on the table in the kitchen, six nodes are placed 
on boards extended outside the windows of the living room. The other twenty nodes are all placed on 
polyvinyl chloride (PVC) stands outside the house. All thirty-four nodes are programmed with TinyOS 
program Spin |[T9l . and a basestation connected to a laptop is used to collect pairwise RSS measurements 
from these nodes. 

Both experiments are performed using the following procedure. Before people start to walk in the 
living room, a calibration is performed with no people (no extrinsic motion) in the experimental area. 
The duration of the calibration period of Experiment 1 is about 47 seconds, and M = 140 measurements 
are recorded for each link; while for Experiment 2, M = 170 measurements are recorded for each link 
during a 57 second calibration period. Compared to L = 1122 directional links, M is much smaller than 
L. Next, a person walks around a marked path in the living room at a constant speed, using a metronome 
and a metered path so that the position of the person at any particular time is known. 

These two through-wall experiments use the same hardware and software, and are performed following 
the same procedure. However, the main difference between these two experiments is the season. Exper- 
iment 1 is performed on a clear winter day, while Experiment 2 is performed on a windy day in late 
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(a) Experiment 1 



(b) Experiment 2 



Fig. 3: Pictures of two experiments. 
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Fig. 4: Experimental layout of Experiment 2. The shade area is covered by tree branches and leaves. 



spring. As shown in Figure 3a there are no leaves on branches and no wind is observed from the video 



of Experiment 1. However, from the video recorded during Experiment 2 (one snapshot is shown in 
Figure [3b]), we observe that wind causes grass, leaves and branches to sway ll20ll . The wind also causes 
the PVC stands supporting the nodes to move. The swaying of leaves and branches and the movement 
of the PVC stands are intrinsic parts of the environment, which cannot be avoided, even when no people 
are present in the environment. Thus, the difference between Experiments 1 and 2 is that Experiment 2 
has more intrinsic motion. 
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IV. Results 



A. Eigen-network results 



As described in Section II-B each of the principal components used to construct the intrinsic subspace is 
an eigenvector of the covariance matrix of the network measurements, and each element in an eigenvector 
is from an individual link, we refer these eigenvectors Uj as "eigen-networks". 

Since the first eigen-network ui = [uu, ni2, • • • , un]^ points in the direction of the maximum variance 
of the calibration measurements yc, we show the first eigen-network ui graphically in Figure [5] We see 
the links with un values higher than 30% of the maximum value are all in the lower right side of the 
house. This is consistent with our observation that the intrinsic motion of the leaves and branches on the 
tree located to the right side of the house causes significant variations in the RSS measured on links likely 
to have RF propagation through the branches and leaves. Note that links with high un values all have at 
least one end point near the tree. In particular, links which are likely to see significant diffraction around 
the bottom-right comer of the house have high un values. The leaves and branches almost touch this 



comer, as seen in Figure 3b Not only do these links measure high RSS variance during the calibration 
period, they do so simultaneously. That is, the fact that these links have high positive un values indicates 
that when one of these links experiences increased RSS variance, the other links also measure increased 
RSS variance. Thus, the first eigen-network ui becomes a spatial signature for intrinsic motion-induced 
RSS variance. When we see this linear combination in y^, we should attribute it to intrinsic, rather than 
extrinsic motion. These observations about the source of RSS variance on links support the idea that 
intrinsic motion in the environment causes increased RSS variance simultaneously on multiple links. 
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Fig. 5: First eigen-network: Links with un > 30% of max^ui/. 
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B. Localization results 

Now, we evaluate VRTI, SubVRT and LSVRT using measurements from Experiments 1 and 2. From 
these three estimators, we obtain reconstructed motion images, and the position of the moving person 
can be estimated by finding the center coordinate of the voxel with maximum value. Specifically, a 
localization estimate is defined as: 



z = Zq where q = arg max Xp 

p 

where Zq is the center coordinate of voxel q, and Xp is the pth element of the estimate x = [xi, ^p]^ 
from Q, ( [IT] ) or ( 14i. Then, the localization error is defined as: eioc = ||z — z\\i^, where z is the actual 
position of the person, and I2 indicates the Euclidean norm. 

The VRTI estimates of Experiment 2 are shown in Figure |6] For clarity, we only show the ac- 
tual/estimated positions when the person walks the last round of the square. We find that due to the impact 
of intrinsic motion, some estimates of VRTI are greatly biased to the right side of the experimental area 
(i.e., five estimates with more than 4.0 m error, as shown in Figure [6]l. However, for SubVRT and LSVRT, 
the impact of intrinsic motion is greatly reduced. As shown in Figure |7] and Figure [8j the estimates from 
SubVRT and LSVRT are more accurate than VRTI. There are no estimate errors larger than 2.0 m. Note 
that for both VRTI and SubVRT, some estimates are outside the house. The algorithms presented do 
not include any prior information of the house map or physical barriers which would prevent certain 
trajectories. Incorporation of prior knowledge of an indoor environment might be used to obtain better 
estimates, but at the expense of requiring more information to deploy the system. 
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Fig. 6: Estimates from VRTI. 
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Fig. 7: Estimates from SubVRT. 
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Fig. 8: Estimates from LSVRT. 



Qualitatively, we compare the localization errors from VRTI, SubVRT and LSVRT for the full data set. 
The comparison between VRTI and SubVRT is shown in Figure |9] and the comparison between VRTI 



and LSVRT is shown in Figure 10 The localization errors from SubVRT are all below 1.8 m. For VRTI, 
there are several estimates with errors above 3.0 m. These large errors are due to the impact of intrinsic 
motion on static link measurements. Specifically, we compare the localization errors during a period with 
strong wind, from sample index 205 to 221, as shown in the inset of Figure [9] During this period, the 
average localization error from VRTI is 3.0 m, while the average error from SubVRT is 0.62 m, a 79% 
improvement, and for LSVRT, it is only 0.50 m, a 83% improvement. 

We also compare the RMSE of the estimates, which is defined as the square root of the average squared 
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localization error over the course of the entire experiment. The RMSEs from the two experiments are 
summarized in Table [I] For Experiment 1, the RMSE from VRTI is 0.70 m, while the RMSE from 
SubVRT is 0.65 m and the RMSE from LSVRT is 0.63 m. Since there are not much intrinsic motion in 
Experiment 1, the improvement in RMSE from SubVRT is 7.0%, and the improvement from LSVRT is 
9.6%. For Experiment 2, the RMSE from VRTI is 1.26 m, while SubVRT and LSVRT are more robust 
to impact of intrinsic motion. The RMSE from SubVRT is 0.74 m, a 41.3% improvement, and the RMSE 
from LSVRT is 0.69 m, a 45.3% improvement. 
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Fig. 9: Estimate errors from VRTI and SubVRT. 
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Fig. 10: Estimate errors from VRTI and LSVRT. 
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RMSE 


RMSE 


Improvement 


RMSE 


Improvement 


Exp. 1 


0.70 


0.65 


7.0% 


0.63 


9.6% 


Exp. 2 


1.26 


0.74 


41.3% 


0.69 


45.3% 



TABLE I: Localization RMSEs from VRTL SubVRT and LSVRT. 



C. Discussion 

The parameters that we use in VRTL SubVRT and LSVRT are listed in Table [III Here, we discuss the 
effects of the number of principal components k on the SubVRT localization results. We also discuss the 
effects of the covariance matrix parameter on the performance of LSVRT. 
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Fig. 11: Localization RMSE vs. principal component number k. 



An important parameter for SubVRT is the number of principal components used to construct the 



intrinsic subspace. As discussed in Section H-B the first k components are used to calculate the projection 
matrix for the intrinsic subspace 11/. If A; = 0, 11/ = 0, then Hi = 112, SubVRT is simplified to VRTI. 



The RMSE of SubVRT using a range of k are shown in Figure 11 Since the first eigen-network ui 
captures the strongest intrinsic signal, when k = 1, the RMSE of Experiment 2 decreases substantially 
from 1.26 m to 0.82 m. Since Experiment 1 has less intrinsic motion, the RMSE decreases from 0.70 m 
when /c = to 0.65 m when /c = 4, a less substantial decrease. We note that as k increases, more and 
more information in the measurement is removed, and the RMSE stops decreasing dramatically, and even 
increases, at certain k. This is because when k becomes very large, the information removed also contains 
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a great amount of signal caused by extrinsic (human) motion. Thus, the performance of SubVRT could be 
degraded if k is chosen to be too large. The parameter A; is a tradeoff between removing intrinsic motion 
impact and keeping useful information from extrinsic motion. For experiments without much intrinsic 
motion, such as Experiment 1, we choose a small k. However, for Experiment 2, with strong impact 
from intrinsic motion, we use a large k. As listed in Table |llj we use = 4 and = 36 for Experiment 1 
and 2, respectively. 

3.5 
3.0 
2.5 



1.5 
1.0 

'^ii) = 10 * 10 10 2 10 ' 10" 10' 10^ 

Fig. 12: Localization RMSE vs. cr^. 




An advantage of LSVRT over SubVRT is that LSVRT does not need to change any parameter due to 
changes in the environment, such as parameter k in SubVRT. Thus, we only investigate parameter 



in LSVRT, which plays the same role of the regularization parameter a in SubVRT. From Figure 12 
we see the RMSE from LSVRT reaches the minimum at 0.63 m, when cj^ = 0.001 and m = 4. Similar 
to functions of a shown in Figure 10 of jTJl, the localization RMSEs from LSVRT are also shallow 
functions of o"^ in the range from 10~^ to 10~^. That is, LSVRT is not very sensitive to this regularization 
parameter in a wide range. 

V. Tracking 



In this section, we apply a Kalman filter to the localization estimates shown in Section IV-B to better 
estimate moving people's positions over time. Then, we compare the tracking results from VRTI with 
those from SubVRT and LSVRT, and show that the Kalman filter tracking results from SubVRT and 
LSVRT are more robust to large localization errors. 
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Psr^mctcr 


Value 


Description 




1 nn 

lUU 


Regularization psrcinietcr 




4 




k 


4,36 


Numbers of principal components in Exp. 1, 2 




0.001 


Variance of human motion 


„2 


2 


Process noise parameter 


„2 


5 


Measurement noise parameter 



TABLE II: Parameters in VRTI, SubVRT, LSVRT and Kalman filter. 



A. Kalman filter 

In the state transition model of the Kalman filter, we include both position {P^, Py) and velocity Vy) 
in the Cartesian coordinate system in the state vector s = [Px-, Py, Vx, Vy]^ , and the state transition model 
is: 

s[t] = Gs[t - 1] + w[i] (21) 

where w = [Q,{),Wx-,Wy\^ is the process noise, and G is: 

10 10 
10 1 
10 
1 

The observation inputs r[t] of the Kalman filter are the localization estimates from VRTI, SubVRT or 
LSVRT at time t, and the observation model is: 



G 



(22) 



r[t]= Hs[t\+w[t] 



(23) 



where v = [w^, Vy^ is the measurement noise, and H is: 



H 



10 
10 



(24) 



In the Kalman filter, Vx and Vy are zero-mean Gaussian with variance cj^, Wx and Wy are zero-mean 
Gaussian with variance cr^ [21]. The parameters cj^ and cr^ of the measurement noise and process noise 
are listed in Table [ill 
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Methods 


VRTI 


SubVRT 


LSVRT 


Results 


RMSE 


RMSE 


Improvement 


RMSE 


Improvement 


Exp. 1 


0.66 


0.57 


13.6% 


0.57 


13.6% 


Exp. 2 


1.21 


0.72 


40.5% 


0.66 


45.5% 



TABLE III: Tracking RMSEs from VRTI, SubVRT and LSVRT. 




Fig. 13: CDFs of tracking errors. 



B. Tracking results 

We use the Kalman filter described above to track the positions of the person. The cumulative 



distribution functions (CDFs) of the tracking errors from Experiment 2 are shown in Figure 13 We 
see that the Kalman filter tracking results from VRTI have many more large errors than SubVRT and 
LSVRT. 97% of the tracking errors from VRTI are less than 3.91 m, while 97% of the tracking errors 
from SubVRT are less than 1.36 m, a 65.2% improvement, and 97% of the errors from LSVRT are less 
than 1.15 m, a 70.6% improvement. We use the 97*^^ percentile of errors to show the robustness of the 
tracking algorithm to large errors, and the CDFs show the tracking results from SubVRT and LSVRT 
are more robust to these large errors. 

We also compare the RMSEs of the tracking results from VRTI, SubVRT and LSVRT, which are 



listed in Table III For Experiment 1, the tracking RMSEs from SubVRT and LSVRT are both 0.57 m, 
a 13.6% improvement compared to the RMSE of 0.66 m from VRTI. For Experiment 2, the tracking 
RMSE from SubVRT is reduced by 40.5% to 0.72 m compared to 1.21 m RMSE from VRTI, and the 
RMSE from LSVRT is reduced by 45.5% to 0.66 m. We note that the tracking RMSEs from VRTI, 
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SubVRT and LSVRT of Experiment 2 are both larger than Experiment 1 due to the impact of intrinsic 
motion. However, for VRTI the tracking RMSE from Experiment 2 has a 83.3% increase compared to 
Experiment 1, while for SubVRT and LSVRT, they only increases 26.3% and 15.8%, respectively. The 
tracking RMSEs from SubVRT and LSVRT are more robust to the impact of intrinsic motion. 

Finally, the Kalman filter tracking results of Experiment 2 from SubVRT and LSVRT are shown in 



Figure 14 For Experiment 2 with significant intrinsic motion, the Kalman filter results using SubVRT 
and LSVRT estimates can still track a person with submeter accuracy. 



(a) 



(b) 
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Fig. 14: Kalman filter tracking results of Experiment 2 from SubVRT (a) and LSVRT (b). 



C. Discussion 

In the Kalman filter, the process noise parameter cj^ should be chosen according to the dynamics of 
the movement. For example, for tracking vehicles, o"^ should be set to a large value. The measurement 
noise parameter depends on how accurate the observation inputs are. Here, we choose cr^ based on 
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the speed of moving people in typical homes, and we test the effect of using different on the tracking 
errors. The tracking RMSEs from SubVRT for Experiments 1 and 2 are shown as functions of in 



Figure 15 If cr^ is too large, the Kalman filter gives very small weights to observation inputs. On the other 
hand, for very small measurement noise parameter, the system dynamic model contributes little to the 
Kalman filter. Thus, the RMSE reaches the minimum when an appropriate balance between observation 



inputs and dynamic model is found. We also note from Figure 15 that for both Experiments, the RMSEs 
are shallow functions of in a wide range from 0.001 to 20. That is, if we give sufficient weights to 
the observation inputs, which are the localization estimates from SubVRT and LSVRT, our Kalman filter 
tracking results are not very sensitive to the measurement noise parameter. 




VI. Related work 

DEL using RE sensor networks has potential applications in surveillance for police and firefighters. 
Different measurements and algorithms have been proposed iH, 0, Q, ||5l. Eor RSS-based DEL, 
there are essentially two types of algorithms: fingerprint-based algorithms and model-based algorithms. 
Like fingerprint-based real-time location service (RTLS) systems, fingerprint-based DEL methods use a 
database of training measurements, and estimate people's locations by comparing the measurements during 
the online phase with the training measurements Q, Q, |[22l . Since a separate training measurement 
dataset is necessary, fingerprint-based DEL needs substantial calibration effort. As the number of people 
to be located increases, the training requirement increases exponentially. Model-based algorithms |5il, 
(m, Il23l provide another approach. A forward model is used to relate measurements with unknown 
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people's positions, and the localization problem can be solved as an inverse problem. An advantage of 
a model-based algorithm is that it does not need such training measurements, however, sufficient link 
measurements are necessary to solve the inverse problem. The proposed subspace decomposition and 

least squares methods have been applied to a model-based DFL method - VRTI, and can significantly 
improve the robustness of position estimates. These methods may also be used in fingerprint-based DFL 
methods, but we leave this as a possible future research topic. 

VII. Conclusion 

In this paper, we propose to use subspace decomposition and least squares estimation to reduce noise 
in RSS variance-based device-free locaHzation and tracking. We discuss how intrinsic motion, such as 
moving leaves, increase measured RSS variance in a way that is "noise" to a DFL system. The signal 
caused by intrinsic motion has a spatial signature, which can be removed by the subspace decomposition 
method. We apply the subspace decomposition method to VRTI, a new estimator we call SubVRT. We also 
propose an LSVRT estimator that directly uses the covariance matrix of the measurement to reduce the 
impact of intrinsic motion. Experimental results show that SubVRT and LSVRT can reduce localization 
RMSE by more than 40%. We further apply a Kalman filter on SubVRT and LSVRT estimates for 
tracking. We find the tracking results from SubVRT and LSVRT are much more robust to large errors. 
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